**FULL TITLE** 

ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** 
**NAMES OF EDITORS** 



Simulations of the Magneto-rotational Instability in 
Core-Collapse Supernovae 



M. Obergaulinger, P. Cerda-Duran, and E. Miiller 

Max- Planck- Institut fur Astrophysik, Postfach 1317, D-85741 Garching, 
Germany 

M.A. Aloy 

Departament d'Astronomia i Astrofisica, Universitat de Valencia, 
Edifici d'Investigacio Jeroni Munoz, C/ Dr. Moliner, 50, E-^6100 
Burjassot (Valencia), Spain 

Abstract. We assess the importance of the magneto-rotational instability in 
core-collapse supernovae by an analysis of the growth rates of unstable modes in 
typical post-collapse systems and by numerical simulations of simplified models. 
The interplay of differential rotation and thermal stratification defines different 
instability regimes which we confirm in our simulations. We investigate the ter- 
mination of the growth of the MRI by parasitic instabilities, establish scaling 
laws characterising the termination amplitude, and study the long-term evolu- 
tion of the saturated turbulent state. 



1. Introduction 

At the end of its hydrostatic evolution, the inner core of a star of more than 
about eight solar masses, consisting of about 1.5 solar masses of iron or oxy- 
gen, neon, and magnesium loses support against its self-gravity and collapses to 
supranuclear density p nnc > 2 x 10 14 gem -3 , releasing a gravitational binding 
energy -E grv ~ 10 53 erg. When this density is reached, a shock wave is launched 
and starts to propagate outward. Still within the core, the shock wave stalls 
and turns into a stationary accretion shock as energy is lost when nuclei are 
dissociated in the post-shock region: the prompt explosion fails. 

Most of the gravitational binding energy liberated during collapse leaves 
the star in form of neutrinos, but a small (~ 1%) fraction is used to unbind 
the envelope and power a supernova (SN) explosion. How a fraction of Eg rv 
is transferred to the shock wave and the surrounding envelope is still not fully 
understood, and forms the central issue of SN theory. 

This is not a problem of energy budget but of energy transfer. The gravi- 
tational binding energy of the envelope (~ 10 50 erg) as well as the kinetic and 
electromagnetic energies of a SN (~ 10 49 and ~ 10 51 erg, respectively) are much 
smaller than E gTy . Consequently, the transfer of a small fraction of E gTV suffices 
to revive the shock wave and to trigger the explosion. However, it is extremely 
difficult to tap efficiently the energy carried by weakly interacting neutrinos. 
The mechanisms proposed to account for shock revival and explosion can be 
grouped into the following categories: 
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The spherical neutrino mechanism Heating of the post-shock matter by in- 
terac tions with t he neutrinos diffusing out of the core revives the shock 
(e.g.. lBethdll990l ). This mechanism works for stars around 8 due to a 
steep density profile facilitating shock propagation ( Kitaura et al.1120061 ). 



Hydrodynamic instabilities Without the restriction to spherical symmetry, 
hydrodynamic instabilities can develop: convection in the proto-neutron 
star (PNS) and behind the shock and the standing accretion shock in- 
stability (SASI) can enhance the efficiency of neutrino heating and the 
transfer of binding energy of the accretin g gas to the surrounding matter 
(Blo ndin et al.ll2003l : iFoglizzo et al.ll2007l 1. They may also account for the 



pron ounced asphericitie s of the explosion and the observed pulsar kicks 
(e.g., IS check et al. 2006). While this mechanism has proven successful in 



detailed 2d simulatio ns of stars up to 15 solar masses employing state-of- 



the-art microphysics (jMarek et al.ll2009l ). its robustness and applicability 



to stars of higher mass is still uncertain. 

Energy transfer by waves Apart from neutrinos, waves may carry energy 
from the central region to the surroundings. Unstable g- modes in the PNS 
excited by accreting matter can lead to the emission of acoustic waves 



trave lling upwards and depositing energy near the shock ([Burrows et al 



2006}) . A similar effect can be cause d by Alfven waves e xcited in the con 



vective region of a magnetised PNS ( Suzuki et al. 20081 ). The importance 



of acoustic waves has been suggested in simulations, whereas the heating 
by Alfven waves has not been explored thoroughly. 

Rotational mechanisms Any rotation of the progenitor will be amplified dur- 
ing the collapse; additionally, differential rotation will be created even for 
an initially rigidly rotating progenitor. If sufficiently rapid, rotation may 
provide sufficient additional energy to unbind the post-shock matter. In 
contrast to the energy of neutrinos, rotatio nal energy can be tapp ed fairly 
easily, e.g., by viscosity or magnetic fields ( Thompson et al.ll2005l ). Addi- 



tionally, a differentially rotating magnetise d fluid can be unstable agains t 
the magneto-rotational instability (MRI) (jBalbus Hawlevl fl99ll . Il998l ). 



leading to field amplification, turbulence and enhanced transport of angu- 
lar momentum. 

None of these mechanisms has been able to provide a robust explanation 
for SNe across the entire mass range. Because of the diversity and complexity of 
the underlying physics, the problem is studied using complementary approaches: 
very detailed simulations including the best treatment of the microphysics that 
can be afforded and, on the other hand, simplified models neglecting most of 
these complications for the benefit of covering a larger parameter space of ini- 
tial conditions. The study of (magneto-)rotational mechanism relies currently 
most ly on simplified mode ls. 



Akivama et al. ( 20031 ) have pointed out that the most basic instability con- 



dition of the MRI, i.e., a negative radial gradient of the rotational profile, f2(r) 
is typically fulfilled in rotating post-collapse cores. They have, based on sim- 
plified spherical modelling, estimated growth times of a few milliseconds and 
saturation fields of the order of 10 15 G. Simplified global simulations (e.g., 
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lObergaulinger et al. 20061 ) have shown that fields of this range (reached only for 
already strongly magnetised progenitors with fields of the order 10 12 G) can have 
significant effects on the dynamics of the explosion, e.g., by extracting rotational 
energy from the inner core or launching jet-like outflows. Additionally, signs for 
the growth of the MRI were found in a few of those simulations. Based only 
on such simplified studies, a thorough assessment of the influence of the MRI is 
still not possible. 



2. Methods and Aim of our Study 

The main obstacle for investigating the MRI by detailed microphysical simula- 
tions is the necessity to resolve the wavelength of the most unstable MRI modes, 
Amri ; which scales with the local Alfven velocity, leading to a prohibitive reso- 
lution requirement: for a realistic initial field of 10 9 G, a region of several 100 
kilometres would have to be covered by a 3d grid with a resolution of a few 
centimetres or metres. Consequently, we investigated the properties of the MRI 
in core-collapse SNe using high- resolution simulations of small representative 



domains in the unstable region (jObergaulinger et al.ll2009l ) 



We studied basic properties of the MRI in the linear regime under typical 
conditions of post-collapse cores differing from the well-studied case of accretion 
discs mainly by the spherical geometry and the potential importance of thermal 
stratification exerting either a stabilising or a destabilising effect on the MRI. 
Since the absolute value of the Brunt- Vaisala frequency, N, of the star can 
exceed f2 strongly, the dynamics of the core may not be dominated by standard 
MRI modes but rather show a more complex interplay between the MRI and 
hydromagnetic convection. 

Even if the instability grows fast enough, it can affect the explosion only if 
a sufficiently large mean total (Maxwell plus Reynolds) stress tensor is reached 
at the termination of the growth and maintained during the turbulent saturated 
phase. This is a long-standing issue in MRI theory, and no final solution has 
been found although simulations have identified a number of factors, physical 
and numerical, affecting the amplitudes. We tried to estimate the saturation 
amplitude by an analysis of the non- linear phase of our models. 



3. Results 

3.1. Instability Analysis 

We analysed the dispersion relation of wave-like perturbations of a background 
in rotational equilibrium endowed with a weak magnetic field with given rota- 
tional profile and Brunt- Vaisala frequency, N, i.e., the oscillation frequency of 
a buoyant fluid element (which is imaginary in a convectively unstable layer). 
Different from accretion discs, the main stabilising agent counteracting gravity 
is not rotation but pressure. Thus, the corresponding models possess a sub- 

Keplerian rotation. 

The instability criteria, given by iBalbusI ( 19951 ). allow for two branches of 



unstable modes: buoyant and Alfven modes, corresponding to unstable buoyant 
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oscillations and unstable Alfven waves, respectively. Both branches are stabilised 
by magnetic tension at very short wavelengths. 

Apart from a stable region of strongly positive gradients of entropy or f2, 
these modes dominate in two different regions of the parameter space. For a 
negative entropy gradient and moderate (positive or negative) values of iV 2 , 
corresponding to a small influence of the thermal stratification, Alfven modes 
dominate. They grow most rapidly in a narrow wavelength region around Amri 
and have a vanishing growth rate at infinite wavelength, leading to stringent res- 
olution requirements in numerical simulations. If — iV 2 > exceeds a threshold, 
buoyant modes set in, while Alfven modes disappear. Similar to convection in 
a non-magnetised fluid, they grow rapidly from infinite to short wavelengths. 

Based on the relative importance of these modes, we can distinguish six 
regions in the space given by the Brunt- Vaisala frequency and the rotational 
profile: a stable region, the classic (accretion-disc) MRI, Rayleigh instability, 
magnetoconvection, a region in which the stabilisation of convection due to rapid 
(rigid) rotation is lifted by magnetic fields, and a mixed regime, respectively. 



3.2. Numerical Models 

We performed a large parameter study of simplified 2d and 3d simulations rep- 
resenting small (a few kilometres wide) sections of the equatorial region of a 
post-collapse core with a cylindrical mesh of a resolution between 0.625 and 20 
metres. We added a weak magnetic field in vertical, z, direction, either uni- 
form or sinusoidally varying with a vanishing net flux, to the gas assumed in 
rotational equilibrium with a prescribed rotational profile and entropy gradient. 

To save computational costs, we neglected neutrino transport entirely and 
replaced the complex nuclear equation of state by a simple analytic approxima- 
tion. Our numerical code is based on a conservative formulatio n of the equations 



of id eal MHD within the constrained-transport framework ([Evans &: Hawlev 



1988) and uses high- order methods for reco nstruction (e.g., the monotonicity- 



preserving schemes bv lSuresh Huvnhlll997l ) and approximate Riemann solvers. 

The simulations confirm all regimes relevant to the MRI in SNe, with nu- 
merical growth rates very well in agreement with the analytic predictions. Under 
conditions of a rapidly rotating core, we expect typical MRI growth times of the 
order of a few milliseconds, i.e., sufficiently fast to affect the explosion develop- 
ing on time scales of several tens or hundreds of milliseconds (see the left panel 
of Fig. 13. 2. 1 for the time evolution of a model with vanishing net flux). 

The initial perturbations grow exponentially at the MRI growth rate given 
by linear analysis, omri- Their growth is mediated by coherent large-scale chan- 
nel modes, i.e., axisymmetric p airs of upflows and dow nflows with alternating 
magnetic polarity. As shown by I Goodman Xul ( 19941 ). the channel modes are 



themselves unstable against secondary, parasitic instabilities of Kelvin-Helmholtz 
or tearing-mode type. Our models show the growth of these parasites feeding 
off the channel modes. Their growth rate is a function of the field strength 
of the channels. The higher the field strength of a channel is, the faster grow 
its parasites, until, at a certain field strength, the amplitude of the parasites 
exceeds that of the channel: the MRI growth terminates. In the left panel of 
Fig. 13. 2,| the growth of these parasites can be identified by the strong increase 
of the z-component of the magnetic field near t = 12 ms. Typically, the para- 
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Figure 1. The time evolution of a model with vanishing net flux of the 
magnetic field (left panel) and its spatial structure in the saturated state 
(right panel). The left panel shows the total magnetic energy (black solid 
line), its z-component (green dotted line), and the Maxwell stress component 
Mvjti> = bvjbtj, (blue dash-dotted line) per unit volume as a function of time 
in milliseconds. The right panel shows the computational box; the radial, 
azimuthal, and vertical directions are indicated by the red (pointing towards 
the observer), green (hidden, pointing to the right), and blue (pointing up) 
arrows. The left and right halves of the figure show a volume rendering of the 
toroidal magnetic field, b<f,, (positive and negative values are shown by blue 
and red colours) and of the vorticity of the flow, respectively. 

sitic mode quenching the MRI in our models is of tearing-mode type. This is 
clearly an artefact of our numerics. Since we evolve the system of ideal MHD 
there should be no resistivity, and thus no tearing-mode instability. However, 
the unavoidable numerical resistivity, albeit only at a small level, permits the 
growth of resistive instabilities. 

Combining the MRI growth rates from linear analysis and those of the 
parasites determined in a series of auxiliary simulations, we establish scaling 
laws for the value of the Maxwell stress at the termination of MRI growth. A 
stronger initial field, 60, corresponds to wider channel flows which can sustain a 
stronger field before they are disrupted by parasites. Consequently, the termina- 
tion amplitude scales with the initial field as (6o) 16//7 . A higher grid resolution, 5, 
corresponds to a smaller value of the numerical resistivity, and thus to a slower 
growth of the parasites, leading to a scaling of the termination amplitude as 
£-8/7 _ si ower rotation leads to wider channels, again decreasing the growth of 
parasites. Hence we find a scaling with Finally, the growth rate of the 

parasites depends on the local sound speed, yielding a scaling with 

After the disruption of the coherent channel flows, turbulence develops (see 
Fig. 13. 24 right panel, for the structure of a model in this state). The flow and the 
magnetic field (shown in the right and left halves of the panel) are characterised 
by a multitude of small-scale filamentary features. Most obvious are magnetic 
flux tubes elongated in azimuthal direction. Depending on the MRI regime 
and the predominance of Alfven (buoyant) modes, efficient transport of angular 
momentum (entropy) leads to rigid rotation (flat entropy profile, and a rotation 
profile where the specific angular momentum is constant). 

The Maxwell stress is typically maintained at a level of the order of its 
termination value throughout saturation, but it is subject to pronounced oscil- 
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lations and even intermittent episodes of exponential growth. These correspond 
to the re-appearance of c oherent channel-like flow s that terminate after a short 
time by parasitic modes ( Sano &: Inutsukall200ll ). Apart from these flows, we 



identify large-scale correlations of the velocity and magnetic field in the turbu- 
lence maintained over several orbital periods, similar to the ones observed by 
iLesur fe Qgilviel (|2008l b 

Our results suggest that the MRI may play a role, for SNe that have rapidly 
rotating progenitors. The field strengths reached in our models are of the order 
of 10 15 G, i.e., sufficiently strong to affect the explosion by redistributing angular 
momentum and extracting energy from rotation. For slower rotation, the inter- 
play of the MRI and convection, i.e., the mixed instability regime may play a 
role in the unstable regions of the PNS and the surrounding neutrino-heated hot 
bubble. Neglecting a large part of the relevant physics and restricting ourselves 
to a rather special geometry, we are, however, far from a final answer, e.g., in 
the form of a MRI model with predictive power for generic post-collapse cores. 
Further modelling is required to include the effects neglected thus far. 
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